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Abstract 

Angular correlations in —)■ A(3872)A+ decays, with A(3872) —)■ 

A —^ vr''“7r“ and J/V' —)■ are used to measure orbital angular momentum 

contributions and to determine the value of the X (3872) meson. The data cor¬ 
respond to an integrated luminosity of 3.0 fb“^ of proton-proton collisions collected 
with the LHCb detector. This determination, for the first time performed without 
assuming a value for the orbital angular momentum, confirms the quantum numbers 
to be The A(3872) is found to decay predominantly through S wave 

and an upper limit of 4% at 95% C.L. is set on the D-wave contribution. 
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The X(3872) state was discovered in 5+’° —)■ X(3872)i^+’°, X(3872) —)■ 7r’''7r J/'0, J/'ip —)■ 
decays by the Belle experiment Q and subsequently confirmed by other experiments 
HQ Its production was also studied at the LHC |^|^. However, the nature of this state 
remains unclear. The X(3872) state is narrow, has a mass very close to the D^D*^ threshold 
and decays to p^J/ip and uJ/tj) final states with comparable branching fractions [^, thus 
violating isospin symmetry. This suggests that the X (3872) particle may not be a simple 
cc state, and exotic states such as D^D*^ molecules [^, tetraquarks or mixtures of 
states have been proposed to explain its composition. The X(3872) quantum numbers, 
such as total angular momentum J, parity P, and charge-conjugation C, impose constraints 
on the theoretical models of this state. The orbital angular momentum L in the X(3872) 
decay may also provide information on its internal structure. 

Observations of the X(3872) —)■ 7J/0 and X(3872) —)■ 7'0(2S') decays 


Ilf 13 imply 


positive O, which requires the total angular momentum of the dipion system (Jtttt) in 
X(3872) —)■ 7r~^7T~ J/'ijj decays to be odd. Dipion mass, is limited by the available 

phase space to be less than 775 MeV, and so > 3 can be ruled out since there are 
no known or predicted mesons with such high spins at such low masses p| In fact, the 


in line 


distribution of is consistent with X(3872) —)■ p^J/'ijj decays [b, 14,15 

with JjTTT = 1, the only plausible value. 

The choices for were narrowed down to two possibilities, l^"*" or 2~~^, by the 
CDF collaboration, via an analysis of the angular correlations in inclusively reconstructed 
X(3872) —)■ and J/'0 —>■ decays, dominated by prompt production in pp 

collisions |16 . Using 1.0 fb~^ of pp collision data collected by LHCb, = 2 was 
ruled out in favor of the 1 assignment, using the angular correlations in the same 
decay chain, with the X(3872) state produced in B~^ —)■ X(3872)X’'' decays [^. Both 
angular analyses assumed that the lowest orbital angular momentum between the X(3872) 
decay products {L^in) dominated the matrix element. Significant contributions from 
Lrain + 2 amplitudes could invalidate the I’*"’'' assignment. Since the phase-space limit on 
is close to the pole (775.3 ± 0.3 MeV [^), the energy release in the X(3872) 
decay, Q = M{J/'ijj7r~^TT~) — M{J/ijj) — is a small fraction of the X(3872) mass, 

making the orbital angular momentum barrier effective]^ However, an exotic component 
in X(3872) could induce contributions from higher orbital angular momentum for models 
in which the size of the X(3872) state is substantially larger than the compact sizes 
of the charmonium states. Therefore, it is important to probe the X(3872) spin-parity 
without any assumptions about L. A determination of the magnitude of contributions 
from Lmin -|- 2 amplitudes for the correct is also of interest, since a substantial value 
would suggest an anomalously large size of the X(3872) state. In this article, we extend 


our previous analysis 17 of five-dimensional angular correlations in —)■ X(3872)X’'', 


X(3872) —)■ p^J/4’, ^ TT+TT”, J/pj —)■ decays to accomplish these goals. The 

integrated luminosity of the data sample has been tripled by adding 8 TeV pp collision 
data collected in 2012. 


^ The inclusion of charge-conjugate states is implied in this article. 
^ We use mass and momentum units in which c = 1. 

^ Dimuon candidates are constrained to the known J/b mass [^. 
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Figure 1: Distribution of AM for —>■ J/'ipK^7r^7r~ candidates. The fit of the A(3872) signal 
is displayed. The solid (blue), dashed (red), and dotted (green) lines represent the total fit, signal 
component, and background component, respectively. 


The LHCb detector is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, described in detail in Ref. (^[^. The X(3872) candidate selection, 
which is based on reconstructing —)• {J/'ijj —t 'k~ candidates using particle 

identification information, transverse momentum (pt) thresholds and requiring separation 
of tracks and the vertex from the primary pp interaction vertex, is improved relative to 
that of Ref. 17 . The signal efficiency is increased by lowering requirements on px for rnuons 


from 0.90 to 0.55 GeV and for hadrons from 0.25 to 0.20 GeV. The background is further 
suppressed without signihcant loss of signal by requiring Q < 250 MeV. The X(3872) 
mass resolution (uam) is improved from about 5.5 MeV to 2.8 MeV by constraining the 
B^ candidate to its known mass and requiring its momentum to point to a pp collision 
vertex in the kinematic fit of its decay. The distribution of AM = J/'ip) — M[J/'ip) 

is shown in Fig. A Crystal Ball function with symmetric tails is used to model 
the signal shape, while the background is assumed to be linear. An unbinned maximum 
likelihood ht yields 1011 ± 38 B^ —)■ X(3872)A+ decays and 1468 ± 44 background 
entries in the 725 < AM < 825 MeV range used in the angular analysis. The signal 
purity is 80% within 2.5 a^M from the signal peak. From studying the A+vr+vr” mass 
distribution, the dominant source of background is found to be B~^ —)■ J/ipKi{1270)~^, 
Ai(1270)+ —j- A+tt+tt” decays. 

Angular correlations in the B~^ decay chain are analyzed using an unbinned maximum- 
likelihood fit to determine the A(3872) quantum numbers and orbital angular momentum 
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in its decay. The probability density function {V) for each hypothesis, Jx, is 
dehned in the hve-dimensional angular space Vt = (cos^x, cos6'p, cos^j/^, A0x,j/p), 

where 6*x, Op and Oj/^ are the helicity angles 21-23 in X(3872), and J/ijj decays, 


respectively, and A0x,p, A0x,j/v> are the angles between the decay planes of the X(3872) 
particle and of its decay products. The quantity V is the normalized product of the 
expected decay matrix element {M.) squared and of the reconstruction efficiency (e), 
VWx) = |M(b!|Jx)Pe(b!)//(Jx), where J(Jx) = / |Xl(b!| Jx)P e(b!)db!. The efficiency 
is averaged over the mass using a simulation (24| - [^ of the X(3872) —)■ p^J/ip, 

—)■ 7r+7r“ decay. The lineshape of the resonance can change slightly depending on the 
X(3872) spin hypothesis. The effect on e(fl) is very small and is neglected. The angular 
correlations are obtained using the helicity formalism (Id), 
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where A are particle helicities, AA^ = A^+ — A^-, and are Wigner functions 

The helicity couplings, Axj^^^x^, are expressed in terms of the LS couplings, B^s, witK the 
help of Clebsch-Gordan coefficients, where L is the orbital angular momentum between 
the and the J/'ip mesons, and S is the sum of their spins. 
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Possible values of L are constrained by parity conservation, Px = Pj/p Pp (—1)^ = (—l)"^- 

1, only the minimal value of the angular momentum, Lmin, 


In the previous analyses 14|16, 
was allowed. Thus, for the preferred = I’*"’*' hypothesis, the D wave was neglected 
allowing only S-wave decays. In this work all L values are allowed in Eq. (§. The 
corresponding B^s amplitudes are listed in Table Values of Jx up to four are analyzed. 
Since the orbital angular momentum in the B^ decay equals Jx, high values are suppressed 
by the angular momentum barrier. In fact, the highest observed spin of any resonance 
produced in B decays is three 29,30 . Since V is insensitive to the overall normalization 


of the Bls couplings and to the phase of the matrix element, the B^s amplitude with 
the lowest L and S is set to the arbitrary reference value (1,0). The set of other possible 
complex Bls amplitudes, which are free parameters in the fit, is denoted as a. 

The function to be minimized is —21nT(Jx,Q!) = — 2 tCj lnP(Oi| Jx, «), 

where T( Jx,«) is the unbinned likelihood, and iV^ata is the number of selected candidates. 
The background is subtracted using the sPlot technique 31 by assigning a weight, Wi, to 


each candidate based on its AM value (see Fig. [^. No correlations between AM and 11 are 
observed. Prompt production of X(3872) in pp collisions gives negligible contribution to 
the selected sample. Statistical fluctuations in the background subtraction are taken into 
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Table 1: Parity-allowed LS couplings in the X(3872) —)• decay. For comparison, we also 

list a subset of these couplings corresponding to the lowest L, used in the previous determinations 
of the X(3872) quantum numbers. 
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JPC 

Bls 

Any L value 

Minimal L value 

0 -+ 

Bn 

Bn 

0 ++ 

Boo, B22 

Boo 

1 -+ 

Bioi Bn, Bi2, Bo2 

Bio, Bin Bi2 

1 ++ 

Boi, B21, B22 

Boi 

2 -+ 

Bill Bi2, Boi, H32 

Bin Bi2 

2 ++ 

Bo 2 -i B20, B21, B22, B42 

Bo2 

3-+ 

Bi2, Boo, Boi, B02, B02 

Bi2 

3++ 

B21, B22, B41, B42 

B21, B22 

4 -+ 

Boi, Bo2, Boi, B02 

Boi, Bo2 

4++ 

B22, B40, B41, B42, Bq 2 

B22 


account in the log-likelihood value via a constant scaling factor, 

The efficiency e(f2) is not determined on an event-by-event basis, since it cancels in the 
likelihood ratio except for the normalization integrals. A large sample of simulated events, 
with uniform angular distributions, passed through a full simulation of the detection and 
the data selection process, is used to carry out the integration, I{Jx) oc x)\‘^ ■, 

where A^mc is the number of reconstructed simulated events. The negative log-likelihood 
is minimized for each Jx value with respect to free B^s couplings, yielding their estimated 
set of values a. Hereinafter, C{Jx) = 'C(Jx,d). 

The hypothesis gives the highest likelihood value. From angular momentum and 
parity conservation, there are two possible values of orbital angular momentum in the 
X(3872) decay for this value, L = 0 or 2. For the S-wave decay, the total spin of 
the and mesons must be S' = 1; thus Bqi is the only possible LS amplitude. For 
the D-wave decay, two values are possible. S' = 1 or 2, corresponding to the amplitudes 
B 21 and B 22 , respectively. The squared magnitudes of both of these D-wave amplitudes 
are consistent with zero, as demonstrated by the ratios |H 2 iP/|i?oiP = 0.002 ± 0.004 
and |i? 22 p/|.BoiP = 0.007 ± 0.008. Overall, the D-wave signihcance is only 0.8 standard 
deviations as obtained by applying Wilks theorem to the ratio of the likelihood values with 
the D-wave amplitudes hoated in the ht and with them hxed to zero. The total D-wave 
fraction depends on the B^s amplitudes, /d = / \A4{Q)i)f dVt/ J |A4(D)g+D|^ dD, where 
A4(D )D is the matrix element restricted to the B 21 and B 22 amplitudes only and A 4 (D) 5 +d 
is the full matrix element. To set an upper limit on /d, we populate the four-dimensional 
space of complex B 21 and B 22 parameters with uniformly distributed points in a large 
region around the B 21 and B 22 ht values (±14 standard deviations in each parameter). For 
each point we determine the likelihood value from the data and an /d value via numerical 
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Figure 2: Likelihood-weighted distribution of the D-wave fraction. The distribution is normalized 
to unity. 


integration of the matrix element squared. The distribution of /d values weighted by the 
likelihood values is shown in Fig. It peaks at 0.4% with a non-Gaussian tail at higher 
values. An upper limit of /d < 4% at 95% C.L. is determined using a Bayesian approach. 

The likelihood ratio t = —21n[£(J^^*)/£(l++)] is used as a test variable to discriminate 
between the I’*"'*' and alternative spin hypotheses considered (Jx*)- The values of t in 
the data (tdata) are positive, favoring the I’*"’'' assignment. They are incompatible with 
the distributions of t observed in experiments simulated under various hypotheses, 
as illustrated in Fig. To quantify these disagreements we calculate the approximate 
signihcance of rejection (p-value) of (fdata — where (t) and a{t) are the 

mean and rms deviations of the t distribution under the hypothesis. In all spin 
conhgurations tested, we exclude the alternative spin hypothesis with a signihcance of 
more than 16 standard deviations. Values of t in data are consistent with those expected in 
the case as shown in Fig. with fractions of simulated experiments with t < fdata 
in the 25%-91% range. Projections of the data and of the ht V onto individual angles 
show good consistency with the assignment as illustrated in Fig. Inconsistency with 
the other assignments is apparent when correlations between various angles are exploited. 
For example, the data projection onto cos 6x is consistent only with the I’*"’*' ht projection 
after requiring | cos0p| > 0.6 (see Fig. |^, while inconsistency with the other quantum 
number assignments is less clear without the cos6p requirement. 
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Figure 3: Distributions of the test statistic t = —2 ln[£(JJ^'*))/£(1++)], for the simulated 
experiments under the hypothesis (blue solid histograms) and under the = I"’"'*' 

hypothesis (red dashed histograms). The values of the test statistics for the data, tdatai 
shown by the solid vertical lines. 


The selection criteria are varied to probe for possible biases from the background 
subtraction and the efficiency corrections. By requiring Q < 0.1 GeV, the background 
level is reduced by more than a factor of two, while losing only 20% of the signal. By 
tightening the requirements on the pt of vr, K and /i candidates, we decrease the signal 
efficiency by around 75% with similar reduction in the background level. In all cases, the 
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significance of the rejection of the disfavored hypotheses is compatible with that expected 
from the simulation. Likewise, the best £t /d values determined for these subsamples of 
data change within the expected statistical fluctuations and remain consistent with the 
upper limit we have set. 



Figure 4: Background-subtracted distributions of all angles for the data (points with error bars) 
and for the 1"’“"’' fit projections (solid histograms). 


In summary, the analysis of the angular correlations in —)■ X[?>'S12)K^, X(3872) —>■ 

TT^'K~ J/tjj, J/'ip —)■ decays, performed for the hrst time without any assumption about 

the orbital angular momentum in the X(3872) decay, conhrms that the eigenvalues of total 
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Figure 5: Background-subtracted distribution of cos^x for candidates with | cos0p| > 0.6 for the 
data (points with error bars) compared to the expected distributions for various X(3872) 
assignments (solid histograms) with the Bis amplitudes obtained by the fit to the data in the 
five-dimensional angular space. The fit displays are normalized to the observed number of the 
signal events in the full angular phase space. 


angular momentum, parity and charge-conjugation of the X(3872) state are I’*'’*'. These 
quantum numbers are consistent with those predicted by the molecular or tetraquark 
models and with the Xci(2^Pi) charmonium state 32 , possibly mixed with a molecule 10 


Other charmonium states are excluded. No significant D-wave fraction is found, with an 
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upper limit of 4% at 95% C.L. The S-wave dominance is expected in the charmonium 
or tetraquark models, in which the X(3872) state has a compact size. An extended size, 
as that predicted by the molecular model, implies more favorable conditions for the D 
wave. However, conclnsive discrimination among models is difhcnlt becanse qnantitative 
predictions are not available. 
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